Collective flow effects on charge balance correlations and local parity-violation 
observables in ^/snn = 2.76 TeV Pb+Pb collisions at the LHC 
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We study the effects of collective flow on charge dependent azimuthal correlations at LHC energies. 
We propose a series of correlations as a signature of the combined effects of azimuthal collective 
flow and local charge conservation and perform an analysis within a statistical freeze-out model. 
We find that present LHC measurements of charge dependent azimuthal correlations are consistent 
with local charge conservation on the kinetic freeze-out surface. In view of experimental searches 
for signatures of the chiral- magnetic effect, we provide an alternative explanation of the charge 
dependence of the observed signal and propose additional measurements to disentangle the effects. 



I. INTRODUCTION 

In relativistic heavy-ion collisions a large number of 
charged particles is produced throughout the dynamical 
evolution of the system. While this is readily evident 
from measurements of the hadronic spectra in the 
final state of the evolution, the spectra provide only 
little insight into how and when the production of 
charged particles occurs. In contrast measurements of 
charged particle correlations contain additional infor- 
mation about the evolution of the system and several 
observables have been proposed to study in particular 
the chemical evolution of the system [l|-(j|. In this 
context it is important to realize, that charged particle 
production is subject to microscopic conservation laws, 
which require local production of charge anti-charge 
pairs in coordinate space. The strong collective expan- 
sion of the system transforms these correlations from 
coordinate space to momentum space, while at the 
same time diffusive interactions with the medium reduce 
the correlations 0, Charge balance correlations 

therefore contain important information on the creation 
mechanisms of charged particles and their subsequent 
evolution. 

In the literature it is usually assumed that charge 
production happens at two different stages of the evolu- 
tion [lHflj. The first wave of production is expected to 
occur during the first fm/c of the collision by initial hard 
scatterings and during the thermalization process to 
form the quark-gluon plasma phase . After a collective 
expansion of the system, a second stage of charge 
production is expected during the transition from the 
deconfined QGP phase to hadronic matter around 5—10 
fm/c. The quark anti-quark pairs produced early in 



the collision arc subject to diffusion over a large period 
of time and one expects the correlation of balancing 
partners to be significantly reduced in the final state. In 
contrast charge anti-charge pairs which are produced at 
times close to the kinetic freeze-out of the system exhibit 
diffusion only over a short period of time, and one 
expects a strong final state correlation. Charge balance 
correlations are therefore sensitive to the production 
time of charges as well as the transport properties of the 
medium. While it is an ongoing effort to disentangle the 
different effects we will study here to what extent 

experimental observations at the LHC are consistent 
with highly localized charge conservation at kinetic 
freeze-out. We will compare results from a statistical 
freeze-out model which takes into account local charge 
conservation at freeze-out @, 0] to present LHC data [8[ 
and provide predictions in cases where observables have 
not yet been measured. 

The observable employed in experimental studies 
of balancing charge correlations is the charge balance 
function 



Bcc(pp\p a ) 



_ N C c(pp\p a ) - N cc (pp\p a ) 

dM/dp a 
dMj dp a 



(1) 
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which aims to identify balancing partner charges on 
a statistical basis [l|, Hj]. Here dM/dp a denotes the 
differential charged particle multiplicity and N C c{p/3\p a ) 
is the number of particle pairs where one particle has 
charge c and momentum p a and the other has charge c 
and momentum pp. The balance function B C c{p a ,pp) 
describes the conditional probability to observe a particle 
with charge c and momentum pp, given the observation 
of a particle with opposite charge c and momentum p a . 
Instead of considering the six-dimensional correlation 
function in Eq. ([T]), previous studies have focussed on 
the integrated correlation functions B(A(j>, <fi) and B(Arj) 
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which quantify the separation of balancing charges in 
azimuthal angle A</> and relative pseudorapidity A?y 
[2-0, 19l4l^|. As a consequence of the anisotropy of the 
system in the transverse plane the balance function 
B(A(j), <fi) is sensitive also to the ang le 6 of the pair with 
respect to the reaction plane 0, [l2|. This behavior 
also constitutes a 'background effect' in experimental 
searches for the chiral magnetic effect [13] and has 
triggered and intensive discussion on the interpretation 
of experimental results @, 0, [U, EMI- 

In this paper we focus on the balance function 
B(A(f>, (f>) and study the dependence on the reaction 
plane angle <fi at LHC energies. We employ a statis- 
tical freeze-out model which takes into account local 
charge conservation at kinetic freeze-out, in order to 
estimate the combined effects of collective flow and local 
charge conservation on charge dependent correlation 
functions. In addition to the charge balance function, 
we also consider mixed harmonic azimuthal correlations, 
(cos(m(4> a — 4>i}) — n{(f)p — ^fc))), which correspond to 
different moments of the balance function with respect 
to the fc-th order reaction plane ^k- In this way we 
can investigate not only the separation width of the 
balancing charge but also its variation with respect to 
the event plane, which is sensitive to the flow properties 
of charged particle pairs. One of the correlations we 
investigate, (cos(</> Q + (j>p — 2^np)), has originally been 
proposed in the search of experimental signatures of the 
chiral magnetic effect (CME) [20(. In measurements of 
the STAR and ALICE experiments HG3, 

it was found 

that this correlation shows a strong charge dependence 
which has been argued to be consistent to the CME 
8, 14]. However the combined effects of local charge 
conservation and elliptic flow have been shown to 
account for the charge dependence of the signal observed 
at RHIC [1,0, [13]. We will present results, which show 
that this statement holds also for the LHC measurement 
and we propose to systematically measure a series of 
correlations to disentangle the different effects. 
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FIG. 1. Charged particle yield dN/dpr and elliptic flow 
V2{pt) from the ALICE collaboration [23. 123I] and the blast 
wave model for six different centrality classes. The param- 
eters of the blast wave model are summarized in the above 
table. 



II. FREEZE-OUT MODEL 

We use the freeze-out model developed in [1,0] in order 
to calculate charge dependent correlations in ^/snn = 
2.76 TeV Pb+Pb collisions at the LHC. The model is 
based on conventional blast wave models, where in ad- 
dition local charge conservation on the kinetic freeze-out 
surface is taken into account [1, 0] ■ Here we use a filled 
ellipse version of the blast wave model, which is a sim- 
ple parameterization of the freeze-out configuration [2ll ] . 
The model parameters of the blast-wave model are the 
kinetic freeze-out temperature Tki n , the flow parameters 
n, (3q, P2, P4, the eccentricity of the freeze-out surface e 
and the chemical freeze-out temperature T c h- The trans- 
verse collective flow velocity is parameterized as 

/3(r, fa) = r"(0i)(Aj + /32cos(20 b ) + &cos(40 6 )) (2) 



The model parameters are entirely determined by single 
particle spectra and are chosen to reproduce LHC results 
of z>2, V4 and pt spectra from the ALICE collaboration 
22l 23 1 . The parameters are summarized in the above 
table and the model is compared to the ALICE data in 
Fig. [T] In addition to the single particle parameters we 
consider the parameter o^, which control the separation 
of balancing charges on the freeze-out surface 7]. We 
will mostly consider the case = 0, where balancing 
partner charges are emitted from the same point on the 
freeze-out surface. This amounts to perfectly local charge 
conservation and gives rise to the strongest possible cor- 
relations. In some cases we vary the parameter to 
study the sensitivity of our results. The implementation 
of local charge conservation is identical to Refs. [1,0] and 
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FIG. 2. The balance function with respect to the reaction 
plane calculated from the freeze-out model. One observes that 
the balance funtion is narrower for in-plane pairs as compared 
to out-of plane pairs. The balance function at intermediate 
angles shows an asymmetry in A(/> towards the reaction plane. 



proceeds as follows: In the first step of calculation we cre- 
ate arrays of particles according to a canonical ensemble 
with fixed volume (100 fm 3 ) and chemical freeze-out tem- 
perature T c h = 148 MeV [23], such that electric charge, 
strangeness and baryon number all sum to zero. We 
then randomly assign a point on the freeze-out surface 
to each array, which determines the collective flow veloc- 
ity. Finally the momentum of each particle is assigned 
according to a thermal distribution with the collective 
velocity of the array. Charge dependent correlations can 
be calculated separately for each array as particles from 
different arrays are assumed to be uncorrelated in our 
model. We note that in reality, charge dependent corre- 
lations are sensitive also to the correlations induced by 
final-state interactions and identical particle interference 
f24j. When incorporated in the model, these have been 
shown to explain distortions of the balance function at 
RHIC However, for our purpose of estimating the 

magnitude of the local charge conservation signal, we ex- 
pect these effects to be less relevant and we will neglect 
them in the following. 



III. CHARGE DEPENDENT AZIMUTHAL 
CORRELATIONS AT LHC ENERGY 

Here we present results for charge dependent azimuthal 
correlations from the blast wave model described in the 
previous section. In order to study the combined effects 
of collective flow and local charge conservation, we first 
investigate the charge balance function in azimuthal an- 
gle and with respect to the reaction plane 
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dM/d(j> 
S(A<p- 



J dpa^-dpp B(p a \pp) 

(<j> a - fo))6(<j> - fa) . (3) 



FIG. 3. Illustration of the combined effects of collective flow 
and local charge conservation. The top panel shows how in- 
plane pairs are more tightly correlated in azimuthal angle as 
compared to out of plane pairs, due to larger collective veloc- 
ities in-plane. To lower panel illustrates that at intermediate 
angles the balancing partner charge is more likely to be found 
towards the in-plane direction. 



Here <f> is the azimuthal angle of the leading particle 
with respect to the reaction plane and A<f> is the relative 
angle of the pair. The results are shown in Fig. [2~J where 
we present the balance function B(A<f>, <p) as a function 
of A(j> for <j) ~ 0° , 45° and 90° and we use the same 
kinematic cuts (\r]\ < 0.8, 0.2GeV/cp T < 5GeV/c) as 
the ALICE collaboration. From Fig. [2]one observes that 
the balance function in-plane (<j> ~ 0°) is significantly 
narrower as compared to the balance function out-of- 
plane (cf> ~ 90°). This is a consequence of the larger 
collective velocities in-plane, which lead to a focussing in 
azimuthal angle as illustrated in the top panel of Fig. [3] 
At intermediate angles (</> ~ 45°) the balance function 
shows an asymmetry towards the in-plane direction. 
This can be attributed to the fact that the balancing 
partner charge is statistically more likely to be found 
in-plane as compared to out-of-plane as illustrated in 
the lower panel of Fig. [31 

The modulation of the balance function in the 
azimuthal angle of the pair <f>, can be quantified in 
terms of harmonic moments of the balance function. In 
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FIG. 4. The centrality dependence of the integrated corre- 
lation <5{cos(Ac/>)) calculated by the freeze-out model. The 
different values of a$ correspond to different separations of 
balancing partners on the freeze-out surface. Black circles are 
the results from the ALICE collaboration. The results from 
HIJING simulations (green dot line), which do not include 
local charge conservation clearly underestimate the observed 
signal. 
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FIG. 5. The centrality dependence of the integrated corre- 
lation 5(cos(raAc/>)} calculated by the freeze-out model with 
o,f, — 0.0. The correlators quantify higher harmonic moments 
of the charge balance function. 



order to establish this relation we introduce the set of 
correlators 



C a p(n,m; k) 



M a M p 



(4) 



and we will only consider the charge dependent part 

5C ee i [C+_ + C.+ - C++ - C__] . (5) 



Here ipk denotes the fc-th order event plane. Since in 
the statistical model the event planes of different orders 
coincide, we will restrict ourselves to k = 2 in the fol- 
lowing. However in the presence of flow fluctuations this 
is no longer the case and it would be interesting to also 
study the behavior in event- by- event simulations. The 
correlators SC(n, m; 2) are related to the charge balance 
function by 



(M 2 5C(n, m; 2)) 

W) 



x cos(mA0 + nq 



(6) 



such that C(n, to; 2) quantifies the n-th harmonic modu- 
lation of the rn-th harmonic moments of balance function. 
The zeroth moment <5C(0, 1;2) = <5(cos(Ac/>)) quantifies 
the overall width of the charge balance function and is 
shown in Fig. The higher (unmodulated) moments 
(5C(0,to;2) = (5(cos(mA</>)) are shown in Fig. [5l In or- 
der to compare our model results to ALICE data, we ac- 
count for the detector efficiency and acceptance by using 
the same kinematic cuts as the ALICE experiment and 
rescaling our data to reproduce the experimental multi- 
plicities and the overall normalization [25[ 

Z = {N+--N++)/N+ + (N-+-N—)/N- . (7) 

One observes that for central and mid-central collisions 
the ALICE results are compatible with perfectly local 
charge conservation on the freeze-out surface. For more 
peripheral collisions we find that local charge conserva- 
tion over-estimates the signal. However if we reduce the 
degree of correlation by taking into account a non-zero 
separation of balancing charges on the freeze-out 
surface, we find that the model results are compatible 
with the ALICE data. We note that final-state inter- 
actions and identical particle interference might further 
reduce the correlation, such that charge conservation is 
more local than suggested by Fig. |U This effect can be 
investigated by measurements of the balance function 
B{4>) as well as measurements of the higher moments 
(5C(0,to;2). For comparison we also show results from 
HIJING simulations, which do not account for local 
charge conservation and clearly underestimate the signal. 

The modulated moments n > generically fea- 
ture three different contributions, which can be written 
as BHEI 



5C(n, m; 2) = v n C m + v$ - 



(8) 



where we introduced the flow coefficients v n and the mo- 
ments of the balance function C m as 



v„ 

Ci-n 



(cos(n0)) , 
<5(cos(toA</>)) 



(9) 
(10) 



The contribution v n C m is a combination of the harmonic 
modulation of the single particle spectra and the charge 
balance width. For instance for the second harmonic mo- 
ment 5C(2,1;2), the contribution V2C1 originates from 
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FIG. 6. The centrality dependence of the integrated corre- 
lation SC(n, 1; 2) calculated by the freeze-out model. The 
correlation 5C(2, 1; 2) has been measured by the ALICE col- 
laboration in the context of local strong parity violation. We 
find that local charge conservation provides an alternative ex- 
planation of the charge dependent signal. 

higher particle multiplicites in-plane as compared to out- 
of-plane 0, 0, HH • The other contributions take the form 

v$ = (cos(mA^)cos(n^)) - v H C m , (11) 
»W = (sin(mA0)sin(n0)) . (12) 

which are associated to the modulation of the charge bal- 
ance function in the azimuthal angle of the leading parti- 
cle. Concerning the second harmonic moment SC(2, 1; 2) 
one finds that the correlator vi. c describes the modula- 
tion of the width of the balance function, whereas V2, s is 
negative and originates from the asymmetry of the bal- 
ance function at intermediate angles [6j, [D, [l5| . These fea- 
tures can clearly be observed from Fig. [2] and similarly 
one can identify the contributions for n — 4,6,... asso- 
ciated to higher order flow harmonics. The blast-wave 
results for the correlation function SC(n, 1;2) are pre- 
sented in Fig. [6] as a function of centrality. The correlator 
SC(2, 1; 2) has originally been proposed as a signature of 
the chiral magnetic effect [2(| and been measured by the 
ALICE collaboration The results in Fig. [6] show that 
local charge conservation already accounts for the entire 
charge dependence of the observed signal. We emphasize 
that this observation is in complete analogy to the situ- 
ation at RHIC [1,0, E3- Fr om Fig. 1 one also observes 
a large signal also for the n — ±4 harmonic components, 
which originates from the quadrangular anisotropy and 
we suggest to use this correlator to verify the origin of the 
correlation. On the other hand correlations with n = ±6 
are highly suppressed because f^l^} is very small. In 



the context of the experimental search for signatures of 
the chiral magnetic effect, it might therefore be interest- 
ing to investigate the higher moments (\n\ > 4), where 
the background contribution from local charge conserva- 
tion is significantly reduced. 

IV. SUMMARY AND DISCUSSION 

In this paper we studied the combined effects of 
collective flow and local charge conservation on charge 
dependent correlation functions at LHC energies. We 
employed the blast wave model from Ref. [a, 0] and 
performed a fit to single particle spectra to extract 
the blast wave parameters. Within this model we 
studied charge balance correlations in the transverse 
plane. Similar to the results at RHIC energies 0, 0, [l2j , 
we find that collective flow leads to a modulation of 
charge balance correlations with respect to the reaction 
plane. This can be quantified in terms of the moments 
SC(n, to; k) of the balance function which we propose to 
measure. When comparing our model results to LHC 
measurements in the context of searches for the chiral 
magnetic effect Q, we find that the charge dependence 
observed in the ALICE data can be fully explained with 
local charge conservation at kinetic freeze-out. We also 
provided predictions for additional correlators to test 
this alternative explanation of the signal. We emphasize 
that the notion of local charge conservation at kinetic 
freeze-out is an interesting observation by itself. This 
is because the correlations of balancing charges in the 
final state are directly related to the chemical evolution 
of the plasma as well as the transport properties of the 
medium. The key challenge is to identify the origin 
of local charge conservation at freeze-out which can 
be explained by small in-medium diffusion of quarks 
produced in the early stage of the evolution as well as 
by a different production mechanism where charged 
particles are predominantly produced at late times. In 
this context it was pointed out that balance functions 
of identified charged particles provide a promising tool 
to disentangle hadronic and partonic contributions 
and one should also take into account the recent results 
of direct photon measurements, which point to strong 
collective behavior early in the collision 27i]. Finally it 
would be interesting to study the effects of local charge 
conservation in event-by-event simulations and possibly 
also within a microscopic approach. First results in this 
direction have been presented in (TH [l9| and we expect 
more studies in the future. 
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